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Abstract 

We theoretically investigate the possibility of realizing a nonlinear 
all-optical diode by using the unique features of quasiperiodic ID pho- 
tonic crystals. The interplay between the intrinsic spatial asymmetry 
in odd-order Thue-Morse lattices and Kerr nonlinearity, combined with 
the unconventional field localization properties of this class of quasiperi- 
odic sequences, gives rise to sharp resonances that can be used to give a 
polarization-insensitive, nonreciprocal propagation with a contrast close 
to unity for low optical intensities. 

1 Introduction and motivations 

An all-optical diode (AOD) is a spatially nonreciprocal device that in the ideal 
case, and for a specific wavelength A, allows the total transmission of light 
along the forward direction (T^~ = 1), and totally inhibits light propagation 
along the backward direction (T^~ = 0), yielding a unitary contrast C = (T^ — 
T^)/(T A h + T^) — 1. AODs are widely considered to be the key components 
for the next generation of all-optical signal processing, in complete analogy 
with electronic diodes which are widely used in computers for the processing 
of electric signals. Replacing relatively slow electrons with photons as carriers 
of information would substantially increase the speed and the bandwidth of 
telecommunication systems, leading to a real revolution of the telecom industry. 
However, controlling photons is a much more difficult task than controlling 
electrons, due to the fact that the latter particles are electrically charged, and 
thus much more suitable for obtaining the nonlinear characteristics typical of 
the electric diode. 

Due to the crucial technological implications, such unidirectional propaga- 
tion ('diode action') has been studied experimentally by several groups, us- 
ing the most diverse schemes and experimental techniques. Photonic crystals 
(PCs) PP, i.e. multilayered structures which possess a definite periodicity of 
the same order of the wavelength of light, are very suitable for this kind of 
applications because of their ability to totally inhibit the propagation of light 
over well-defined frequency regions, known as photonic bandgaps (PBGs, see 
Refs. [U H]). Especially noticeable in this direction are the early efforts of 
Scalora et al. [3], who in 1994 proposed the use of ID PCs, with a gradation 



1 



that makes the structure spatially asymmetric, allowing nonreciprocity due to 
different forward/backward nonlinear shifts of the band edges in the transmis- 
sion spectrum. More recently, in 1999 Gallo and Assanto [I] proposed the use 
of a waveguide with an asymmetrically placed defect, which breaks the for- 
ward/backward symmetry, in presence of quadratic nonlinearities. Most of the 
above schemes, however, suffer from some serious drawbacks which make them 
not suitable for commercial and large-scale applications. Relatively large phys- 
ical sizes are often needed; the balance between figures of merit and optical 
intensities is usually inadequate, and in some cases cumbersome structural de- 
signs are necessary to provide structural asymmetry. For instance, in the widely 
used optical isolators, which allow a strong unidirectional propagation and are 
based on the Faraday rotation effect [5], the use of two polarizers is necessary, 
because the effect is strongly polarization-dependent, which makes the device 
very efficient but a few centimeter long. 

In this work, we discuss a different type of AOD design based on ID nonlin- 
ear multilayered structures arranged in a quasiperiodic or aperiodic fashion [6] . 
Since the first report by Schechtman et al. [7] of metallic Al-Mn alloys showing a 
quasiperiodic crystal structure, and the pioneering experimental works of Merlin 
et al. on Fibonacci and Thue-Morse (T-M) GaAs- AlAs superlattices [5] , a large 
number of studies have been devoted to ID, 2D and 3D quasicrystals, see [1] 
and references therein. More specifically, we examine the diode action in T-M 
lattices, which exhibit isolated, sharp and high quality resonances inside the so- 
called pseudobandgaps, i.e. bandgap regions which are not due to translational 
symmetry. Thanks to the intrinsic asymmetry of odd-order T-M sequences, and 
due to the unique field localization properties of T-M sequences, we have found 
that it is possible to achieve nonrcciprocal nonlinear shifts of these 'pseudores- 
onances' by using small optical intensities, which combined with the small size 
of the structures (of the order of a few microns), would make them feasible for 
advanced applications and for integration on chip. 

2 Linear transmission of T-M photonic crystals 

A quasiperiodic chain is defined by a recursive, deterministic generation rule 
So = {A}, A — > fi(A, B) and B — > f%{A, B), where So is the initial string and 
fi^(A, B) are two arbitrary strings of symbols A and B. A and B indicate the 
two kinds of layer of which the structure is made of, with refractive indices ua 
and n b respectively. Generation Sj is a string which is constructed by applying 
the above recursive rule j times. Several basic quasiperiodic chains have been 
extensively investigated in the literature, such as the Fibonacci, Cantor, Double- 
Period and Rudin-Shapiro sequences, see also [5] and references therein. The 
quasiperiodic structure that we consider here, and one of the most important 
in general, is the T-M sequence, which was introduced by Thue in a 1906 study 
[9]. It is based on the rule T = {A}, A -> {AB}, B {BA}. The size of the 
string Tj for the T-M sequence is given by 2 J , so that for real- life applications 
one usually considers j < 10. The first T-M sequences are: T = {AB}, T — 
{ABBA}, T 3 = {ABBABAAB}, etc. The ratio L 3 between number of blocks of 
type A and number of blocks of type B in Tj is always equal to unity, irrespective 
of the generation number j. It is important to note that the T-M sequence Tj 
is symmetric with respect to string inversion operation T (which is performed 
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by reversing the order of the string elements) if j is an even integer, while it 
is antisymmetric (i.e. it is invariant under X plus the substitution A — ► B 
and B — ► A) if j is an odd integer. In the present work we consider only 
antisymmetric sequences, in that asymmetry, together with nonlincarity, is an 
essential requisite for the nonreciprocal behavior which is at heart of the optical 
diode action. 

Figs. l(a,b,c,d) show the linear transmission spectra for plane waves (as 
calculated with the transfer matrix method, TMM, see Ref. [2]) for the T- 
M photonic crystals ^3,5,7,9 respectively, in the wavelength range 0.5-1.1 //m. 
Normal incidence is assumed, i.e. vanishing incidence angle 6 — 0, so that 
TE- and TM-polarization spectra are degenerate. At this stage, forward and 
backward spectra also coincide (T& = TC, VA), since the linear transfer ma- 
trix M is unimodular [2J, and, despite the strong asymmetry of the considered 
lattices, no optical diode action is possible in the linear regime. Crystal pa- 
rameters are: refractive indices of the layers ua — 1-55 (which corresponds to 
the polydiacetylene 9-BCMU organic material, see Refs. |T0l nB = 2.3 
(which corresponds to TiC"2 material), An = ns — Ua — 0.75, layer thicknesses 
d-A = Ao/4/n^i ~ 112.9 nm and ds — A /4/ris ~ 76.1 nm, where Ao — 0.7 
/j,m (indicated with a dashed red line in Fig. 1). Total physical lengths of the 
crystals are approximately 0.756 fj,m for T3 [8 layers], 3.02 /jm for % [32 layers], 
12.09 fj,m for T 7 [128 layers] and 48.38 fim for T 9 [512 layers]. 

Several things can be deduced from the panel of Fig. 1. First of all, it is easy 
to see that, when increasing the generation number j, each peak of transmission 
shows evidence of self-similarity. This is a very well-known feature which is 
discussed in many works [llj . Secondly, most of the peaks shown in Fig. 1 are 
peaks of perfect transmission (PPTs, for which T = 1), regardless of the value 
of j, and they tend to accumulate around the reference frequency Ao, at which 
the quarter wavelength condition is satisfied. We have found that, amongst the 
basic quasiperiodic crystals, this unexpected property seems to be a unique and 
remarkable feature of T-M sequences. Another important feature shown in Fig. 
1 is the gradual (for j > 7) formation of pseudobandgaps, i.e. frequency regions 
of zero transmission (not due to periodicity) where the crystal behaves like a 
perfect mirror. In Figs. 1(c) and 1(d) one particular pseudobandgap under 
scrutiny is indicated with a green shaded region, while other pseudobandgaps 
(the ones containing other PRs) present in the spectra are indicated with blue 
areas. This gap contains a sharp resonance [pseudoresonance, PR, indicated 
with blue arrows in Fig. 1(c) and 1(d)], analogous to a defect state, that in 
presence of nonlinearity can be used to obtain nonreciprocal propagation. This 
is located at A res ~ 0.8643 /im. The generation of the above sharp resonances 
can be attributed to the specific localization phenomena in T-M sequences [T2J . 
Both extended and critical states exist in T-M sequences, which give rise to two 
distinct types of bandgaps. One is the traditional bandgap that exists in periodic 
crystals, the other one shows fractal (self-similar) splitting into smaller bands. 
The above types of gaps are distinguished by their different behavior when 
increasing the generation number [12j . Traditional gaps change their depths 
only, while fractal gaps split, creating the sharp peaks of perfect transmission 
which are impossible to obtain with periodic crystals or resonators. 

Quality factor (Q-factor) is a measure of the average energy in the resonant 
peak over the energy radiated per cycle [T5]. It is defined by Q = \ res /AX, where 
A res is the central resonant wavelength and AA is the resonance full-width at half 
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maximum (FWHM). In Fig. 2(a) the evolution of Q-factor for the PR indicated 
above is shown as a function of generation number. The Q-factors for the PRj 
resonance are: Q(PR 6 ) = 262.06, Q(PR 7 ) = 864.8, Q(PR&) = 1736.53 and 
Q(PRq) = 1235.4. The resonance appears (together with the pseudobandgap) 
at j = 6 [Fig. 2(b)], and it has a well-defined bell-shape until j = 9 [Fig. 
2(e)], at which the fractal self-replication of the spectrum occurs. At this point 
the PR broadens due to the birth of satellite peaks. Remarkably, the central 
wavelength of the PR (where it exists) does not depend on j. For the chosen 
index contrast, it thus makes sense to talk about this PR only for j = 6, 7, 8. 
However, we need to restrict ourselves to odd-order T-M sequences, the only 
ones that do not exhibit spatial symmetry. For this reason we select Tj- as our 
representative structure, which has the best PR specifics for the selected value 
of An. In order for the pulse to clearly distinguish the resonance, its duration t p 
must be ^ 2tt/Auj, where Aw is the resonance FWHM. In the case of PR7 we 
have Aw ~ 2.5 THz, so that a quasi-CW pulse is defined by t p > 2.5 psec. This, 
together with material relaxation time, will limit the device switching time. 

Before embarking in the nonlinear analysis of the structures, it is useful to 
understand how small random fluctuations in the angle of incidence of light from 
perfect normal incidence (9 = 0) affect the transmission spectra for TE- and 
TM-polarizations. Fig. 3 shows the ^-dependence of TE- and TM-polarization 
spectra for T7. Position of the PR7 peak is also indicated. For 9 = (normal 
incidence), the two output spectra must coincide. The most important feature 
that can be deduced from Fig. 3 is the presence of minima in the peak positions 
in correspondence of the 'degeneracy' point 9 = 0, a desirable feature for what 
concerns device applications. This leads to a solid structural stability of the 
PR to small random fluctuations in 9, which will therefore be polarization- 
independent to a high degree of approximation. 

3 Nonlinear transmission and optical diode ac- 
tion 

We now proceed to analyze the nonlinear behavior of the PR7 peak shown in 
Fig. 2(c). In order to do this, we use the nonlinear TMM method outlined 
in Ref. [3J. According to this method, for each value of the input wavelength 
A, we must guess the output forward amplitude, while the backward output 
amplitude vanishes. After propagating back to the input layer, one obtains a 
value for the total input intensity, which is proportional to the sum squared 
of the input forward and backward fields. We repeat this process until the 
input intensity matches the desired incident intensity. Figs. 4(a,b) display a 
sketch of the different nonlinear frequency shifts of PR7 for three different in- 
put intensity conditions (/ = 0, black dotted lines; / = 10 MW/cm 2 , solid 
red lines; and / = 50 MW/cm 2 , solid blue lines), for both forward [Fig. 4(a), 
T?] and backward [Fig. 4(b), TT] incidence. The nonlinear Kerr refractive 
indices for TiC>2 and 9-BCMU used in our calculations are roughly estimated 
to be n 2 (Ti0 2 ) ~ 10" 14 cm 2 /W and ra 2 (9-BCMU) ~ 2.5 • 10~ n cm 2 /W. These 
are measured values 1141 110] , obtained by using degenerate and non-degenerate 
pump-probe experiments. Results published in [lOj show that linear and non- 
linear absorptions can be minimized considerably in 9-BCMU inside certain 
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spectral regions, which are therefore suitable for device operation. Thus here 
we neglect these two absorptions for simplicity. Xhe relaxation time T re i of 9- 
BCMU is estimated to be around 1.6 psec. The switching time of the device is 
fast and equal to max{r p , r re i}. 

Let us take, for instance, an input pump intensity equal to / = I p = 10 
MW/cm 2 , corresponding to the curves displayed with red solid lines in Figs. 
4(a,b). Note the strong difference in transmission between forward and back- 
ward incidence. Moreover, there is a reduction of the transmission peaks with 
increasing intensity. By tuning the pump wavelength near the maximum of 
the forward incidence peak [X p ~ 0.8659 /im, see green vertical dashed line in 
Fig. 4(a)], one has substantial transmission for forward incidence of pump light 
(T+ = 0.905). However, for the same wavelength, the backward transmission 
is located on the tail [see green vertical dashed line in Fig. 4(b)], which ensures 
a very low transmission for the backward direction (TT = 0.06). For the par- 
ticular case shown in Figs. 4(a,b), this gives an efficient diode action with an 
AOD contrast equal to C = 0.875. The resonant nature of the PR enhances the 
nonlinear effects to a point that the AOD action is obtained by using modest 
intensities. 

4 Pulse propagation 

Finally, in order to confirm the results obtained above for realistic quasi-CW 
pulses, which have a small but finite spectral width well inside the resonance, 
we accurately solve Maxwell's equations in presence of an inhomogeneous per- 
mittivity e(z). Let us write the ID Maxwell equations in the conventional 
form dt[eE] + cd z B = 0, d t B + cd z E = 0, where E(z,t) and B{z,t) arc the 
scalar electric and magnetic fields respectively, c is the speed of light in vac- 
uum. e(z,t) = e (z) + e 2 I is the dielectric permittivity, where eo(z) 1 / 2 = n(z) 
is the refractive index (real and > 1), and e 2 is the Kerr-nonlinear dielec- 
tric constant. / = \E(z,t)\ 2 is the electric field intensity. We then decouple 
the fast spatiotemporal oscillations due to the carrier wavenumber kg and fre- 
quency uj with the expansion: E(z,t) = [A(z,t)e ik ° z ~ iuJot + c.c.]/2, B(z,t) = 
[G(z,t)e lk ° z ~ iuJot + c.c.]/2. A and G are the complex envelopes of the electric 
and magnetic fields respectively, and fco = ujq/c. After defining new dimension- 
less variables £ = z/zq and r = t/to, where zq and to are respectively arbitrary 
spatial and temporal scales, we choose z = ct and t = 2n/u . We also rescale 
the fields as {^4, G} — > Aq{A,G}, where Aq is a convenient reference intensity. 
This leads to the dimensionless equations 



which are suitable for numerical computations, and do not rely on any slow- 
varying amplitude approximation. An important observation is that the non- 
linear dielectric function e is given in the new variables by c(£,t) = £q(£) + 
(3/4)e2^4ol^(£i t )| 2 j a non-trivial expression that satisfies simultaneously the 
terms rotating with the first and the third harmonics. Eqs. (jT]|2]) are solved by 
first initializing A as a function of £, at r = 0, and then posing G(£, r = 0) = 
±7j(£)A(£,t = 0), where the plus (minus) sign refers to a forward (backward) 
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propagation condition, which is used to test the forward/backward asymmetry 
in our T-M structure. 

Results of simulations are shown in Fig. 4(c,d). Fig. 4(c) shows the nonlinear 
transmission of a forward supergaussian pulse with a temporal width of r p = 2 
psec and a spatial width £ p = 600 /im (3> than the structure's length, which is 
12.09 /zm), and input wavelength A p = 0.8659 /im, corresponding to the location 
of the green dashed line in Figs. 4(a,b). The chosen supergaussian profile well 
approximates the incidence of a finite CW wave, and input intensity was taken 
to be / = 10 MW/cm 2 . According to the previous discussions, this allows almost 
unitary transmission (T+ w 0.9) if the spectral width of the pulse is narrower 
than that of the PR. Launching the same pulse in the backward direction we have 
a dramatically different behavior: Fig. 4(d) shows the nonlinear transmission in 
this case, where light is almost entirely reflected (T~ ss 0.05), confirming both 
qualitatively and quantitatively the optical diode action and the plane-wave 
predictions discussed above. 

5 Conclusions 

In conclusion, we have studied the linear and nonlinear transmission properties 
of a micron-size T-M quasiperiodic PC near a PR, i.e. a sharp resonance located 
in the pscudobandgap. Sharp resonances are due to the unique field localization 
properties of T-M sequences. We have shown that the strong asymmetry of odd- 
order Thue-Morse lattices, combined with a Kerr nonlinearity, gives rise to a 
highly nonreciprocal transmission which is the major feature of an AOD. This 
effect has also been confirmed by numerically integrating Maxwell's equations 
for realistic optical pulses. Our design allows an unprecedented reduction in 
size of the device at relatively low operational optical intensities, a consequence 
of the intrinsic antisymmetry of the structures considered, and the resonant, 
strongly localized nature of the transmission. 

The author would like to thank Victor Grigoriev for useful comments. This 
work is supported by the UK Engineering and Physical Sciences Research Coun- 
cil (EPSRC). 
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Fig. 1 Linear transmission spectra for plane waves at normal incidence 
(6 = 0) for asymmetric odd-order T-M photonic crystals (a) T$, (b) T$, (c) 
T 7 , (d) Tg. Reference wavelength Ao = 0.7 yum is indicated with a red dashed 
line. The pseudobandgap and the PR of interest for this paper are indicated 
with a green area and a blue arrow respectively Other large pseudobandgaps 
which contain other PRs are indicated with blue areas. Crystal parameters are: 
riA = 1.55, n B = 2.3, d A = \ /4/riA ^ 112.9 nm and d B = A /4/n B ~ 76.1 
nm. 

Fig. 2 (a) Q-factors of PR at A res — 0.8643 /im as a function of generation 
number for T-M PCs. (b,c,d,e) Transmission spectra of the X res — 0.8643 /im 
PR for j = 6, 7, 8, 9 respectively. Crystal parameters are the same as in Fig. 1. 

Fig. 3 ^-dependence of TM- and TE-polarization spectra (left and right 
parts of the plot respectively) for the linear T-M quasicrystal T-j. Large blue 
regions correspond to pseudobandgaps. Crystal parameters are the same as in 
Fig. 1. 

Fig. 4 (a,b) Sketch of nonlinear shifts for forward [(a)] and backward 
[(b)] incidence, for three different input intensities. (c,d) Results of pulsed 
simulations, showing the asymmetric transmission properties of the nonlinear 
device for forward [(c)] and backward [(d)] incidence. In both cases, input 
intensity is I = 10 MW/cm 2 . The input pulse has a supergaussian profile, 
A(£,t = 0) = v / ?exp[— (C/Cp)" 1 ], with m = 4 and £ p = 600 /im, which corre- 
sponds to a pulse duration of t p = 2 psec. The position of the T-M crystal is 
indicated by a vertical bar. Crystal parameters are the same as in Fig. 1. 
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